A novel approach in the WIMP quest: 
Cross— Correlation of Gamma— Ray Anisotropics and Cosmic Shear 
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Both cosmic shear and cosmological gamma-ray emission stem from the presence of Dark Matter 
(DM) in the Universe: DM structures are responsible for the bending of light in the weak lens- 
ing regime and those same objects can emit gamma-rays, either because they host astrophysical 
sources (active galactic nuclei or star-forming galaxies) or directly by DM annihilations (or decays, 
depending on the properties of the DM particle). Such gamma-rays should therefore exhibit strong 
correlation with the cosmic shear signal. In this Letter, we compute the cross-correlation angular 
power spectrum of cosmic shear and gamma-rays produced by the annihilation/decay of Weakly 
Interacting Massive Particle (WIMP) DM, as well as from astrophysical sources. We show that this 
observable provides novel information on the composition of the Extra-galactic Gamma-ray Back- 
ground (EGB), since the amplitude and shape of the cross-correlation signal strongly depends on 
which class of source is responsible for the gamma-ray emission. If the DM contribution to the EGB 
is significant (at least in a definite energy range), although compatible with current observational 
bounds, its strong correlation with the cosmic shear makes such signal potentially detectable by 
combining Fermi-LAT data with forthcoming galaxy surveys, like Dark Energy Survey and Euclid. 
At the same time, the same signal would demonstrate that the weak lensing observables are indeed 
due to particle DM matter and not to possible modifications of General Relativity. 

PACS numbers: 95.35.+d,95.30.Sf,98.62.Sb,98.80.-k,95.85.Pw 



Introduction. Weak gravitational lensing refers to the 
small distortions of images of distant galaxies, produced 
by the distribution of matter located between galax- 
ies and the observer Contrary to strong lensing, 
where distortion is caused by a small number of mas- 
sive objects, weak lensing studies aim at reproducing 
the statistical properties of the density field acting as a 
lens, as well as investigating the geometrical properties 
of the Universe. A distorted image can be described by 
the so-called distortion matrix, normally parametrized 
in terms of the convergence k (controlling modifications 
in the size of the image) and the shear 7 (accounting for 
shape distortions). Whilst the former is a direct estima- 
tor of matter density fluctuations along the line of sight, 
the latter is easier to measure, through correlations in 
the observed source ellipticities. In the flat-sky approx- 
imation, the two generate identical angular power spec- 
tra: we shall thus focus on the shear as an estimator of 
the convergence, and, from now on, we indicate with k 
the observables related to cosmic shear. 

The auto-correlation between the gravitational shear 
in two different directions can provide information on 
the clustering of the large scale structures responsible 
for the lensing effect. The technique has already been 
used in where data of the COSMOS galaxy sur- 

vey allow for a measurement of the two-point corre- 
lation function (2PCF) of the shear at angular scales 
between 0.1 and 10 arcmin. Future surveys, like Pan- 
STARRS, Dark Energy Survey (DES) Q and Euclid 
jToj , due to their larger coverage and improved sensitiv- 



ities, will be able to reconstruct two-dimensional shear 
maps, from which one can extract the auto-correlation 
angular power spectrum (PS). 

The same Dark Matter (DM) structures that act as 
lenses can themselves emit light at various wavelengths, 
including the 7-ray range. While 7-rays can be pro- 
duced by astrophysical sources hosted by DM halos (i.e., 
star- forming galaxies (SFG), or active galactic nuclei 
(AGN)), DM itself may be a source of 7-rays, through 
its self annihilation or decay, depending on the proper- 
ties of the DM particle. Those 7-rays emitted by DM 
should therefore have the potential to exhibit strong 
correlation with the gravitational lensing signal. 

In this Letter we propose to study the cross- 
correlation of gravitational shear with the Extra- 
galactic Gamma-ray Background (EGB), i.e., the resid- 
ual radiation contributed by the cumulative emission 
of unresolved 7-ray sources, as a novel and potentially 
relevant channel of DM investigation. 

The most recent measurement of the EGB was per- 
formed by the Fermi-LAT telescope in [ll[, covering 
a range between 200 MeV and 100 GeV: the emission 
is obtained by subtracting the contribution of resolved 
sources (both point-like and extended) and the Galac- 
tic foreground (due to cosmic rays interaction with the 
interstellar medium) from the whole Fermi-LAT data. 
Unresolved astrophysical sources like blazars (l2- 14|. 



SFGs [15|, [16( or radio galaxies [17|, [18| contribute to 
the EGB but the exact amount of their contribution is 
still unknown. 7-rays produced by DM annihilation or 
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decay can also contribute to EGB |19| - |23| . However, the 
fact that the EGB energy spectrum is compatible with 
a power-law, without any evident spectral feature, sug- 
gests that DM cannot play a leading role in the whole 
energy range 2J, |25[ . In the angular anisotropics of the 
EGB emission, the DM also plays a subdominant role: 
indeed, a detection of a significant auto-correlation an- 
gular PS has been recently reported [26| (for multipoles 
I > 100, which is the range of interest for our analysis 
since the contamination of the Galactic foreground can 
be neglected), but the features of such a signal (in par- 
ticular its independence on multipole and energy) seem 
to indicate an interpretation in terms of blazars [27l . [28| . 

Both cosmic shear and 7-ray emission depend on the 
large scale structure of the Universe, either because this 
is what generates the lensing effect or because those 
same structures can produce 7-rays, hosting astrophys- 
ical sources, or directly from DM annihilation/decay. A 
certain level of cross-correlation between cosmic shear 
and 7-ray emission is therefore expected. The key point 
of our analysis is to understand whether the shear/7- 
rays cross-correlation is within reach of future galaxy 
surveys, as DES and Euclid, and under which circum- 
stances such signal can be proficiently used to disentan- 
gle a true DM signal from the other astrophysical 7- 
rays sources. It is the first time, to our knowledge, that 
such analysis is performed. If successful, this approach 
could provide direct evidence that what is measured by 
weak lensing surveys is indeed due to DM and is not, 
for instance, a manifestation of alternative theories of 
gravity. We will show in the following that in the case 
of a significant DM contribution to the EGB (at least 
in a definite energy range), although compatible with 
current observational bounds, the strong correlation of 
the 7-rays emission with cosmic shear makes the DM 
signal potentially detectable by combining Fermi-LAT 
measurements with DES and Euclid. A strategy could 
then be envisaged to use the cross-correlation results 
to extract information on the DM particle properties 
(namely its mass and annihilation cross section or de- 
cay lifetime Q). 

On more general grounds, this approach can also pro- 
vide valuable information on the clustering of the mat- 
ter field and on the properties of galaxies, complemen- 
tary to what can be deduced from the study of the shear 
or 7-ray auto-correlations alone. Both the normaliza- 
tion and the shape of the cross-correlation angular PS 
will depend on which class of sources is considered as re- 
sponsible for the 7-ray emission and, as a consequence, 
any detection of cross-power can potentially provide in- 
formation on the composition of the EGB. 

Theoretical modeling. The source intensity along a 
given direction n can be written as: 
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FIG. 1. Left: EGB emission as a function of observed energy 
for the four extragalactic components described in the text. 
Data are from Right: 7-ray angular PS at E > 1 GeV 
for the same models of the left panel. The observed angular 
PS is summarized by the black band [2^ ]. 



where xi z ) is the radial comoving distance, g is the 
density field of the source, and W is the window func- 
tion (which does not depend on n). We then define a 
normalized version W = (g) W, so (I g ) = J d\ W(x)- 
Expanding the intensity fluctuations of two source pop- 
ulations i and j in spherical harmonics, one can compute 
(see Appendix) the cross-correlation angular PS (here 
in the dimensionless form): 
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where f g = [g(x\m, z)/g(z) — 1] (f g is its Fourier trans- 
form) and the Limber approximation (fc = ijx) is as_ 
sumed to hold. We consider the sources to be character- 
ized by a parameter m (typically the mass), and g(x\m) 
is the density field of an object associated to m, while 
g(z) = (g(n,z)). Pij can be computed following the so- 
called halo-model approach. The two-point correlation 
is given by the sum of two components, the 1-halo and 
2-halo terms, i.e. P^ 



plh 

ij 



p2h 

ij 



[H-ii 



dn 



P^{k)= j dm ■^f:(k\m)f j (k\m) 



P?f(k) 



dmi - — &i(mi)/*(fc|mi) 
ami 



dm? 



(hi 

dm? 



bj{m-i)jj(k\m2) 



(3) 



P lin (k) ,(4) 



where dn/dm is the number density distribution of 
sources, P lln is the linear matter PS, and bi{m) is the 
linear bias between the object i and matter. Note that 
the average of (g) is given by: 

f dn f 

g(z) = (g(n, z)) = dm ~^ l \ d 3 xg(x\m,z), (5) 



I g {n) = I dxg(x,n)W(x) , 



(1) 



which implies that at small k (where / ~ 
J d 3 x g{x\m)/g) the terms in the square-brackets in 
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FIG. 2. Left: Window functions vs. redshift. For 7-ray 
sources we consider the flux above 1 GeV normalized to the 
total EGB intensity measured by Fermi-LAT. Right: Three- 
dimensional PS of cross-correlation shear/7-rays at z = 0. 



Eq. are ~ 1 (except in the case of a significant bias). 
The 2-halo term is thus normalized to the standard lin- 
ear matter PS at small k, which motivates the normal- 
ization of the window function introduced above. 

We aim at cross-correlating the shear signal (source i 
in Eqs. (HQ) with 7-rays emitted by DM, SFGs, and 
blazars (source j in Eqs. j2][4j). ^ or wna ^ concerns 
weak lensing, W takes the form (see, e.g., Q): 
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and dN/dx represents the redshift distribution of the 
sources, normalized to unity area (such that we can 
take (Ii) = 1 in Eq. ©). For Euclid dN/dz = 
A £ z 2 e"' z / Zo ' L5 , where z = z m /\A with z rn = 0.9 
being the median redshift of the survey and Ae is 
fixed by the normalization J dzdN/dz — 1. For DES, 
dN/dz = A D (z a + z°^)/(z b + c), with a, b, and c pro- 
vided in Table 1 of [32j, and Ad fixed by the normal- 
ization. Since gravitational lensing is sourced by the 
potential wells of the large scale structure, whose Pois- 
son equation relates to the matter distribution p, we 
have g K {x) = p(x), and f K (x) is given by the density 
contrast 5(x). For the bias in Eq. (|4]) we use the esti- 
mates in [33l | . We adopt the halo mass function dn/dm 
of 34 1, the halo concentration from (35|, and a NFW 
halo density profile [36j . 

For the case of 7-rays from decaying DM we again 
have fd(x) = S(x) (we assume p m ~ pdm)- The win- 
dow function is now given by: 
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where m x and are the mass and decay lifetime of the 
DM particle, J d = /~ dE dNd( f E 1+z)) e -r(E(i+z),z) with 

dNd/dE(E) being the number of 7-ray photons emitted 
per decay in (E E + dE), and r being the optical depth 
for absorption [37[. Note that the factor QdmPc comes 
from the normalization of W, since in this case (g^) = p- 



The DM annihilation signal scales with p 2 , thus we 
have f a = u(k\m) given by the Fourier transform of 
p 2 (x\m) / (p 2 ) . In the literature, equations are often 
written in terms of the so-called clumping factor: 
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where (a a v) is the velocity-averaged annihilation rate 
(which we assume to be the same in all halos) and J a = 
/~ dE^(E(l+z)) e -T(.W+*),*) with dN a /dE{E) be- 
ing the number of 7-ray photons emitted per annihi- 
lation in the energy range (E,E + dE). In the an- 
nihilating DM case, the predictions for both the win- 
dow function and the PS heavily depend on the (un- 
known) clustering at small masses (i.e., on the min- 
imum halo mass, concentration below approximately 
1O 6 M , and on the amount of substructures). We con- 
sider m m j n = 10 _6 Mq (typical free-streaming mass for 
WIMPs) and include unresolved subhalos following the 
scheme described in [38| with parameters tuned as in 
the HIGH scenario of Sec. 3.3 in [23| (within our halo 
model, it induces only moderate boost factor ~ 2). 

The formalism sketched in Eqs. (jT] [5j) can be used 
also for a population of astrophysical sources, by replac- 
ing the mass with the source luminosity C as the char- 
acterizing parameter. This leads to the replacement of 
dm dn/dm with d£§, where $ is the 7-ray luminosity 
function (GLF). For the range of multipoles of interest 
[l < 10 3 ) both blazars and SFG can be approximated as 
point sources and we have gs(£, x — x') = £5 3 (x — x'), 
which leads to: 

(k,z)= f maxi) dC$(C,z) 7 ^ T v(k\m(£)) 
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with v{k\m) being the Fourier transform of p(x\m)/p 
and (gs) = J d£<& £. In Eqs. © a relation between 
the source luminosity £ and the host-halo mass m is re- 
quired. We compute the source bias bs through the halo 
bias by means of bs{£,z) = b n (Mh(£),z), for which 
we need again a relation between host-halo mass and 
source luminosity. On the other hand, since at low red- 
shift and in the mass-range of interest bh ~ 1, the 2-halo 
term is only very mildly dependent on the description 
of M n (£). For a power-law spectrum with index a, the 
window function is: 
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FIG. 3. Cross-correlation of the cosmic shear with the 7-ray 
emission for the four benchmark models described in the 
text (with a 7-ray threshold expected for Fermi-LAT after 5 
years of exposure and a distribution of background sources 
for lensing as in Euclid) . Each contribution is normalized by 
multiplying Eq. @ by (Ij) / (Iegb) to make them additive. 
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FIG. 4. Left: Cross-correlation between cosmic shear and 
7-ray emission, for different classes of 7-ray emitters. DES 
is taken as the reference galaxy survey. Error bars are es- 
timated for the total signal (in black). Right: Same as in 
the left panel but for annihilating DM, with Euclid as the 
reference galaxy survey. 



where Eq = 100 MeV and As is a factor that depends on 
which specific luminosity is chosen as the characterizing 
parameter (as we will describe below). 

The GLF of blazars is computed following the model 
described i n 1331 with the AGN X-ray luminosity func- 
tion from |40| and with the numerical value of pa- 
rameters derived in [iU by fitting Fermi-LAT data on 
EGB diffuse emission and anisotropies. The spectrum is 
taken to be a power-law with a = 2.2, and C is the 7-ray 
luminosity at 100 MeV (which leads to A s = (l + z)~ a ). 
We assume that no blazars fainter than the luminos- 
ity cutoff £ m in = 10 42 erg/s can exist at any redshift, 
while £ max (z) is the maximum luminosity above which 
a blazar can be resolved (for 5-yr Fermi-LAT it is com- 
puted taking F max = 2 • 10 -9 cuT^s -1 for E > 100 
MeV which is the value assumed in Figs. (J2H4J) ) . The 
relation between halo-mass and blazar luminosity can 
be described through M h = lO 11 - 3 M (£/lO 44 - 7 crg/s) 1 - 7 
following [4l| where the blazar 7-ray luminosity is linked 
to the mass of the associated supermassive black hole, 
which is in turn related to the halo mass. The descrip- 
tion of Mh (£) suffers from sizable uncertainties (see Ap- 
pendix) which propagate to the prediction of the 1-halo 
term. However, as we will see later, the blazar contri- 
bution is largely subdominant, thus such uncertainties 
do not affect our conclusions. 

For the GLF of SFGs, we follow results from the 
Fermi-LAT Collaboration [42[, which are based on the 
infrared (IR) luminosity function derived in [43j], and 
the rescaling relation between 7-ray and IR luminosity 
obtained analyzing resolved SFGs [42j . The spectrum is 
assumed to be a power-law with a = 2.7, similarly to the 
Milky Way case, and C is the 7-ray luminosity between 
0.1 and 100 GeV (which leads to A s = (a-2)/(l + z) 2 ). 
The relation Mh(C) could, in principle, be computed 
from the relation between 7-ray luminosity and star 
formation rate (SFR) [42|, the Schmidt-Kennicutt law 
(connecting SFR and gas density), and the ratio of gas 
to total galactic mass. This leads to different relations 



for each different subpopulation of SFGs (e.g., ellipticals 
are much brighter than spirals of the same mass) ; on the 
other hand we do not have 7-ray data to compute the 
specific GLF of the sub-populations, thus we have to de- 
rive an effective averaged relation. Assuming a power- 
law scaling M h = _4-10 12 M Q (£/10 39 erg/s) 8 and a max- 
imum galactic mass of M max = 10 14 M Q , we can find A 
and B using, e.g., the Milky Way data (M h ~ 10 12 M Q 
and £ ~ 10 39 crg/s) and requiring that the mass associ- 
ated to the maximum luminosity ~ 10 43 erg/s (this can 
be computed from the maximum measured IR luminos- 
ity rescaled to 7-ray frequency) not to exceed M max . 
We found A ~ 1 and B ~ 0.5. This is just a simple 
model, and we estimated the impact of the associated 
uncertainty on the 1-halo term (by varying A and B 
within reasonable ranges) to be 0(1) (see Appendix). 

Results. For the sake of clearness we focus on a bench- 
mark annihilating (decaying) DM scenario, where the 
WIMP has a mass of 100 GeV (200 GeV), annihilation 
(decay) rate of (a a v) = 8 • 10- 26 cm 3 /s (r d = 3 ■ 10 26 s) 
and dominant final state bb. The characteristics of the 
DM particle are chosen to saturate (at least in one par- 
ticular energy range) the EGB emission, without violat- 
ing the experimental constraints. In particular, we note 
that, although we take DM to be a significant compo- 
nent of the EGB at E > 1 GeV in Fig. it is basically 
impossible to obtain an evidence for DM from the an- 
gular PS of 7-rays alone because the latter is dominated 
by the blazar contribution. 

In Fig. [2] we show the ingredients of Eq. ([2} for the 
computation of the shear/7-ray cross-correlation angu- 
lar PS: the window function for the cosmic shear sig- 
nal nicely overlaps with the DM window function, both 
for annihilating and decaying DM, while this happens 
only at intermediate redshifts for the SFG window func- 
tion and only at high redshifts for the case of blazars. 
This suggests that a tomographic approach could be a 
powerful strategy to further disentangle different con- 
tributions in the angular PS (this will be pursued in 
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a future work [44| ) . The shear signal is stronger for 
large masses. The same is true also for the 7-ray signal 
from DM and this fact gives a large 1-halo contribu- 
tion which dominates starting from k < 1/i/Mpc in 
Fig. [3b. Galaxies have masses < 1O 14 M , thus they 
correlate with the shear signal of lower-mass halos and 
the 1-halo contribution becomes important at slightly 
smaller scale k > 1/i/Mpc. Since the bulk of unre- 
solved blazars in 5-yr Fermi-LAT will be hosted in rel- 
atively small halos at large redshift, the 1-halo term of 
the blazar/shcar PS is suppressed. Thus, an important 
result is that, since both the shear and DM- induced 
7-ray signals are stronger for large halos, their cross- 
correlation is more effective with respect to the case of 
astrophysical sources. This, together with the sizable 
overlapping of the DM 7-ray and shear window func- 
tions at low redshift, leads to the expectation of a siz- 
able DM signal in the angular PS, which is indeed what 
we find in Fig. [3] For £ < 100 the 2-halo term dominates 
for all the sources, thus the relative size is roughly given 
by the relative contribution in the total EGB emission. 
At £ > 100, the 1-halo term starts to be important in 
the DM case which grows more rapidly than the astro- 
physical sources. At £ > 10 3 , the 1-halo term takes over 
also in the SFG spectrum which is brought again close 
to the DM curve. Blazars are largely subdominant in 
the whole range of multipoles. 

The observational forecasts for the cross-correlation 
between DES or Euclid and Fermi-LAT are shown in 
Fig.H]for the benchmark models considered in this work 
(for error estimates see e.g. 13 ill , and observational per- 
formances are taken from [Illfllll). The fi gure shows 
that, provided the DM is a significant component of the 
total 7-ray EGB as in our assumptions, a DM signal can 
be disentangled in the angular PS at £ < 10 3 . 

Conclusions. In this Letter, we discussed the cross- 
correlation angular power spectrum of weak lensing cos- 
mic shear and 7-rays produced by WIMP annihila- 
tions/decays and astrophysical sources. We showed that 
this method can provide novel information on the com- 
position of the EGB. Since the shear signal is stronger 
for structures of larger masses and most of the 7-ray 
emission from decaying and (in many models of) anni- 
hilating DM is also produced in large mass halos, their 
cross-correlation is typically stronger than the case of 
astrophysical sources (which are associated to galactic- 
mass halos). Forthcoming surveys like DES and Euclid 
can thus potentially provide evidence for WIMPs. 
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General Formalism 

Expanding the intensity fluctuations of Eq. ([I]) in spherical harmonics, SI g (n) = I g (n) — (I g ) = 
{Ig)^i m CL£ m Yi m {ri), and using orthonormal relation for Yi m , one obtains the expression for the (dimensionless) 
aim coefficients: 

aim = jpj J dnSI g (n)Y; m (n) = -ry J dnd X fg(x,r)W{x)Y e * m (n) , (9) 

where f g = g/(g) — f. Then, through the spatial Fourier transformation of f g and making use of the Rayleigh 
formula, we get: 



If dk - 

<H m = j^ J dnd X j^f g (x,ky k - r W(x)Y; m (n) 



J2 if h'(k X )Y i l m ,(k)Y em ,(n) 



I'm' 



= ±_J M d x ^f g (x,k) 

d f r dk ~ 

= — J d x w( x ) J f g ( x ,k)Mk x )Y; m (k) . (io) 

The angular PS is then Cf = (J2 m a lm a lm*)' wnere hj = 1,2 label the signals (e.g., 7-rays and cosmic-shear), 
and Cp 12 ^ = Cf (which can be easily understood from ag = {— l) ,n a|_ m ) 0- Using Eq. (|TU1) . it takes the form: 

cf } = ^-y J d x Wi( x ) j d x ' WM J ^1 J ^1 (/«(x,*)4(x'^)>i!Wif(* , r')nn l (fe(fc') 

dx Wi( X ) ! d X ' Wjix') I dk P l0 (k,x,x')k{kr)3i'(kr')Y lm {k)Yl, m ,{k) 



2 



dX Wi( X ) / d X ' dkk 2 W ] {x')Pi J (k,x,x')jt(.kr)j t {kr') 



r %W i (x)W j ( X )P ij (k=i/x,x)- 



where in the second step we introduced the definition of through (f Qi {x^)fg 3 = (27r) 3 <5 3 (/c — 
X, x')> an d in the last step wc assumed the Limber approximation to hold for such PS. 
So far, we just re-derived the general expression for the cross-correlation angular PS in Eq. (|2|). The next step is 
to compute the form for the specific PS Pij of interest. The PS is the Fourier transform of the 2PCF in real space, 
^ 2 \x — y) = (f gi (x)f gj (y)). To compute it we follow [29| and assume the density field can be written as the sum 
of independent seeds: 

f(x) = f(rrii, x — Xi) = / dm / d 3 x' <5 3 (cc' — Xi)S(m — mj)/(m, x — x') , (11) 
% J J i 

where i labels the seeds, S is the Dirac-delta function, and we take the mass m to be the parameter which 
characterizes the seeds. The seed density is given by dn/dm = (Y]. 8 3 (x — Xi)8(m — mi)), where, here, () denotes 
the ensemble average over of all possible seed distributions. 
Using Eq. (TTTJ), the 2PCF (fi(x)f 2 (y)) is given by: 

£ (2) (a;,y) = J dmi dmi d 3 xi d i x^(^^ i & A {x\ — Xi)S(mi— m,) ^ 5 3 (x2 — Xj)S(rri2— w.j))/i(mi, x — xi) ji{m-2., y— x 2 ) . (12) 



1 Note that if there are k components contributing to a signal /fc = h/Iy }g thc fraction of contribution from each component. 

y (e.g., blazars and DM for 7-rays), one has I y = Ylk ^k an d 
Sly/Iy = Y] k SI k /I y , which leads to af2 = J2 k fk a im> where 
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The correlation of a seed with mass mi at position x\ with a different seed with mass m 2 at position x 2 is provided 
by the seed 2PCF Q (mi,TO2, x\ : x?), and it is not difficult to see that: 

8 z (x-l - Xi)5(mi - mi) 8 z (x 2 - x 3 )S(m 2 - m,)> = 4— -r— [1 + d 2) ( m i> m 2, sbi,»2)1 + -r— -<5 3 (cci - £c 2 )<5(mi - m 2 ) 
* — ' * — ' ami am2 ami 

» i 

(13) 

which leads to: 



/d 
dmd 3 x 1 —f 1 (x - x 1 ,m)f 2 (y -Xi,m) 



dm 1 dm 2 d 3 x 1 d 3 x 2 -^—-^— fi(x - Xi,mi)f 2 (y - x 2 ,m 2 )^ 2 \mi,m 2 ,Xi,x 2 ) . (14) 
ami dm 2 

Remembering that P(k) is the Fourier transform of ^ 2 \x,y) and writing fi in terms of their Fourier transforms 
/i, one obtains: 

Pij{k) = J dm ^f?(k\m) fj(k\m) + J dm 1 dm 2 ^^-^-J*{k\m l )f j {k\m 2 )P {s \k,mi 1 m 2 ), (15) 

where the PS of the seed distribution P^ is the Fourier transform of ^ S 2 K 

Our reference case obviously concerns mass density fluctuations. Therefore Q is the (homogeneous and isotropic) 
linear correlation function of matter (note that it has to be the linear one since in Eq. (|11[) we wrote the density 
field as a linear superposition of seeds). 

For the objects considered in the following (DM halos and astrophysical sources), we will assume (as usu- 
ally done) that their 2PCF ^ can be related to the linear correlation function of matter by means of 

(2) (2) 

Qj (Xi, Xj\mi, m 2 ) w bi(m\) bj(m 2 ) £ l;n (\xi — Xj\) where bi(m) is the linear bias between the object i and matter. 
Thus using P^ (fc, mi, m 2 ) = bi (mi) bj{m 2 ) P hn (k), we finally arrive to Eqs. ([3]) and (j4}. 

Three-dimensional power spectra 

In this section, we explicitly write down and plot the auto- and cross-correlation three-dimensional power spectra, 
which can be derived from Eqs. (|3]) and ((4]) by using the functions gi introduced in the main text. 

Auto- correlation 

For the cases of weak lensing shear and 7-rays from decaying DM, the density field of the source is given by the 
density contrast S(x), i.e., fd{x) = S(x) and f K (x) = 5(x). Therefore the sum of the one and two-halo terms in 
Eqs. (|«|]) and (fj| provides the non- linear matter PS P ss : 



dm^v(k\m) 2 , Pim 



dn. 

dm — — bi l (m)v(k\m) 
dm 



P hn (k) , (16) 



with v(k\m) being the Fourier transform of p(x\m)/p. In order to validate our method, we can thus compare our 
halo-model-based predictions for the matter PS with latest results from high- resolution A-body simulations (46| . 
This is shown in Fig. [5^,. A very good agreement (in particular at low redshift which is the most relevant epoch 
for our purposes) is achieved. 

Approximating astrophysical sources as point sources, we have gs(C,x — x') = £5 3 (x — x'), and Eqs. © and 
(j4|) become: 



P^-ys (k, z) = dC $(£, z) [£^j 4 , P 2h sls (k, z) = £ m ^ Z) dC $(£, z) b s (jC, z) ^ P" n (fc, z) . 

(17) 

Since the term in the square-brackets docs not depend on k (because of the point-source approximation) and Eq. ([5]) 
becomes (gs) = J d£Q £, the 2-halo term is just a rescaled version of P hn with the rescaling factor due to the bias 




FIG. 5. Left: Comparison between the non-linear 3-dimensional matter PS obtained with the halo-model considered in 
this paper and the revised halofit results from high-resolution iV-body simulations of [46| at different redshift. Central: 
Auto-correlation 3-dimensional 7-ray PS of blazars and SFGs at z = 0. For blazars the 1-halo term strongly dominates. 
Right: Auto-correlation 3-dimensional PS of annihilating DM in the three scenarios described in the text at z — 0. 



(so possibly varying with z). The 1-halo term is constant in k, namely it is a "Poisson- noise" term. The results for 
the blazar and SFG benchmark models described in the main text are shown in Fig. 

The signal of annihilating DM scales with p 2 , so to derive the associated angular PS, we can again use the 
formalism outlined in the previous section by just squaring the r.h.s. in Eq. (with / in the r.h.s. still given by 
p). The correlation function would then lead to a term proportional to: 

S 3 (x 1 -x l )5(m 1 -m i ) ^ S 3 (z 2 - x j ) 5(m 2 - rrij ) y 6 3 (x 3 -x k )5(m 3 -m k ) ^ 5 3 (x 4 -xi)5(m 4 -mi)) , (18) 

i j k I 

which can be in turn expanded into four-, three-, and two-point seed correlation functions. However, since halos 
(described by p(x\m)) are mutually exclusive, an annihilation can only occur if the two particles are within the 
same halo. This implies that all the terms other than the 1- and 2-halo contributions (which are the same appearing 
in the r.h.s. of Eq. [TB")) vanish. Eq. (TTJ} does not change, with f a now being p(x\m) 2 / (p 2 ) . The auto-correlation 
PS is thus given by: 



P, 



1/1 



(k) 



dm 



dn ( u(k\m) 



dm 



A 2 



P; 



21, 



(k) 



dn u(k\m) 
dm d^ bh{m) ^- 



P lin (fc) 



(19) 



with u(k\m) being the Fourier transform of p 2 (x\m)jp 2 . Our results are only very mildly dependent on the upper 
mass cutoff (which we set to be m max = 10 18 Mq), while crucial ingredients are the minimum halo mass m min and the 
DM clustering at very small masses. In Fig. [5k, we show three examples: two models having m mm = 1O _6 M0 (which 
is the typical WIMP free-streaming mass) with/ without substructures (the latter described as mentioned in the 
main text) and m m j n = 10 7 AjTq (which is the minimum halo mass currently inferred from dynamical measurements) 
without substructures. 



Cross-correlation 



For the cross-correlation between shear and decaying DM, the 3-dimensional PS P Kld is simply the non-linear 
matter PS P ss . 

The case of cross-correlation of 7-rays from astrophysical sources with decaying DM is analogous to the case of 
7-rays from astrophysical sources with weak lensing shear signal, both given by: 



7S<5 



(k,z) 



«« r 

d£<S>(£,z)——i{k\m{£)) 

c min {z) {9s) 

""^ d£^(£,z)bs(£,z) 7 ^- 



dn 

dm — — bh(m) v(k\m) 
dm 



(20) 

P lin (k,z). (21) 




k[hMpc"'j k[hMpc"'j k[hMpc"'j 

FIG. 6. Left: 3-dimensional cross-correlation PS between shear and 7-rays from blazars at z — 0. The two models for m(C) 
are taken from Ref. [S[ (here C is the 7-ray luminosity at 100 MeV). Central: Same of left panel but for SFGs (with L now 
being the 7-ray luminosity between 100 MeV and 100 GeV). Right: 3-dimensional cross-correlation PS between shear and 
annihilating DM in the three scenarios described in the text at z = 0. 



Results are plotted in Fig. [6] for the benchmark models described in the main text and highlighting the impact 
of m(C) by showing different choices for such relation. We assume the maximum mass for a galaxy to be M max = 
10 14 Af Q . In the case of blazars, the 1-halo term strongly depends on m(C). On the other hand, even with a 
dramatic increase of its level (model B) with respect to our benchmark case (model A), the contribution of blazars 
in the cross-correlation angular PS of Fig. [3] remains subdominant. The dependence of the SFG-shear PS on the 
m(C) relation is milder. 

The cross-correlation power spectra involving the annihilating DM case can be computed from the term 



(^2 S 3 (x 1 - Xi)S(mi - m,) ^ 5 3 (x 2 - x j )5{m 2 - rrij) ^ S 3 (x 3 - x k )S(m 3 - m k )) 



(22) 



an d ag ain noting that only terms in the r.h.s. of Eq. (|T3|) survive , leading again to Eq. (|14[) (see also Appendix B 
in |31|). This leads to: 



dn u(k\m) 
am — — v(k\m) — — — 
am A z 

dm — — bh(m)v(k\m) 
dm 



dn u(k\m) 
dm d^ bh{m) ^ 



P Iin (fc) 



(23) 



and the 3-dimensional PS of cross-correlation of annihilating DM with shear and decaying DM are again analogous. 
Results are shown in Fig. [5] for the same three annihilating DM models of Fig. [5] 

For completeness, although not directly related to the computations shown in this work, we provide the PS in 
the case of cross-correlation of 7-rays between astrophysical sources and annihilating DM: 



" ( 'W,*)A* ( * M£)) 



p2h 

7S7. 



(M 



A 2 
C 



C ma *(z) 

d£$(£,z)b s (£,z) 
c mtn (z) \9s) 



dn u(k\m) 
dm d^ hh{m) ^^ 



P Un (k,z) 



(24) 



